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Abstract 



^ , A bound-state field theory approach to muonic hydrogen is set up using a variant of the Furry 

C"*^ ' representation in which the lowest-order Hamiltonian describes a muon in the presence of a point 

Coulomb field, but the origin of the binding field is taken to be three charged quarks in the proton 
which are modeled as Dirac particles that move freely within a spherical well. Bound-state field 



Oh. 

^ ' theory techniques are used to evaluate one- and two-photon effects. Particular attention is paid to 



two-photon exchange diagrams, which include the effect of proton polarizability. In addition the 



^ , modification of the electromagnetic self energy of the proton by the electric field of the muon is 

t^^ ' examined. Finally, the model is used to carry out a calculation of the static electric polarizability 

o : 

^^ ' of the proton. 
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I. INTRODUCTION 

One of the simplest ways to model the proton is as three very light quarks confined in 
a spherical well. Choosing the radius of the well to be 1.2 fm leads to moderately good 
agreement with experiment for its electromagnetic properties, such as the charge radius, 
magnetic moment, and static electric and magnetic polarizabilities. This model, a simplified 
version of the MIT bag model [l|, will be referred to in the following as the static- well model. 
It allows an alternative approach to the calculation of the electromagnetic properties of the 
proton, generally treated with methods quite different in character, that uses the methods 
of conventional bound-state QED. The latter theory is characterized by wave functions that 
satisfy the Dirac equation in an external field along with electron propagators defined in 
terms of the same field. When the external field is that of a point Coulomb source, a 
modification of the interaction representation introduced by Furry J2J allows a systematic 
Feynman diagram treatment of radiative corrections. This approach can also be applied to 
many-electron systems, and a Feynman diagram treatment of electron-electron interactions 
is also possible. As will be explained below, the present paper is patterned on a calculation 
of these interactions in heliumlike ions involving two-photon exchange 3|. 

The approach we will use in this paper was applied some time ago [4| to the computation 
of the electromagnetic self energy of the proton and neutron. In that work, both the effect 
of exchange of a photon between quarks along with the electromagnetic self energy of the 
quarks were evaluated and found to sum to 0.53 MeV for the proton and —0.28 MeV for the 
neutron for the case of nearly zero-mass quarks. The fact that the proton is lighter than the 
neutron remains explained by the fact that the down quark is heavier than the up quark, 
but it is of note that the electromagnetic correction to the mass splitting, —0.8 MeV, is the 
same order-of-magnitude as the neutron-proton mass difference, 1.2 MeV. 

The proton can be studied with electron-scattering experiments, which have a long his- 
tory of providing information about its properties, in particular the root-mean-square (rms) 
radius, r^. The proton size has recently received considerable attention because of unex- 
pected results for the 2si/2 — 2p3/2 transition energy of muonic hydrogen p]. The issue of 
determining rp from scattering data can be problematic, as extrapolating the slope of the 



bi- 



Dirac form factor to g = involves a number of assumptions [6|. An alternative approach 
is to determine the proton size by doing precise measurements of atomic transitions that are 



sensitive to the effect of the size. The 2010 CODATA result [SJ] in fact uses this procedure 
with hydrogen and deuterium, where the experiment and theory are so accurate that the 
proton size can be inferred with an accuracy comparable to that available from scattering 
experiments as of 2010. 

Because of its smaller size, muonic hydrogen has long been recognized as a system whose 
spectrum could be used to determine a much more accurate rms radius of the proton than 
that obtained from hydrogen and deuterium, but the associated experimental obstacles have 
only recently been overcome. While indeed much more accurate, the result of Ref. (5| for 
the proton size, 

Tp = 0.841 84(67) fm, (1) 

is significantly smaller than the CODATA result, 

Tp = 0.8768(69) fm. (2) 

This discrepancy is referred to as the muonic hydrogen puzzle. 

One possible explanation of the puzzle involves the electromagnetic structure of the pro- 
ton, and the largest theoretical uncertainty comes from an effect called proton polarizability. 
This is generally evaluated by relating the energy shift to forward virtual photon-proton scat- 
tering. The amplitude describing this scattering, T^'^(z/, g^), can then be related to proton 
form factors through dispersion relations. A recent paper that covers all contributions to 







muonic hydrogen with particular attention to proton polarizability is Ref. [7|; in the conclu- 
sion, we compare our results to results quoted in that paper. A number of issues involving 
convergence of the dispersion theory integrals and the need for experimental data compli- 
cate that approach. The purpose of the present paper is to provide an alternative analysis 
patterned after bound-state field theory calculations in atomic physics. This will be done by 
using the static-well model of the proton together with standard bound-state QED. As we 
will show, there is a natural way of setting up a consistent QED calculation for hydrogen and 
muonic hydrogen, with the proton treated as a bound state of three quarks interacting with 
an electron or a muon, that requires no scattering information for its predictions; rather it 
depends only on the radius of the well. 

Regarding the proton as three relativistic particles confined to a small volume is closely 
analogous to treating three electrons in highly- charged ions, where the electrons for large 
nuclear charge Z are quite relativistic and the ion has a size of 1/Z Bohr radius. This 



problem has recently been addressed with techniques similar to those used for heliumlike 
ions mentioned above Isl , and have been shown to provide an accurate description of these 

fin n 

ions [9|, [lOj , the spectra of which have been measured with high accuracy [11 1 . 



In these calculations almost all of the important physics is described by Feynman dia- 
grams with one or two photons. The same turns out to hold for the present calculation, 
though in this paper, while we will show all relevant diagrams, we concentrate our atten- 
tion on two effects dependent on proton structure, the polarizability of the proton and the 
screening of the proton electromagnetic self energy. 

Our model of the proton is extremely simple, but there are three reasons we have chosen 
it. The first is that proton structure effects are generally very small, with even the largest, 
the effect of its finite size, accounting for about 2 percent of the transition energy in muonic 
hydrogen. Thus even a crude determination of a proton structure effect will have a small 
relative theoretical error. The second is that while the proton polarizability correction has 
been evaluated with other methods, a contribution we term the proton Lamb shift has not, 
and the results presented here may stimulate more sophisticated calculations. The final 
reason is that mentioned above, to explore a method of calculating the effect of proton 
structure on atomic energy levels that does not require the use of dispersion theory. 

We will in the following consider the effect of proton structure on both electronic and 
muonic hydrogen. Because our formalism does not include recoil, we will present results in 
terms of the electron mass rrie and the muon mass nip, even though reduced-mass effects on 
the latter are about 10 percent. When we give a general formula, we refer to a lepton with 
mass mi. The state of the lepton, in practice either 2s or 2p3/2, will be denoted v, and the 
index m will be used for sums over intermediate leptonic states: for the corresponding case 
of quarks, we use g to denote a ground-state quark and reserve n for sums over intermediate 
quark states. 

The plan of the paper is as follows. We begin in section II with a quantum mechanical 
(QM) treatment of the shift in energy levels arising from the perturbation of replacing the 
potential of a point proton Coulomb field with that of a general distribution of charge p{r). 
This perturbation theory is evaluated through second order. In section III we turn to a 
quantum-field-theoretic approach to the problem in the context of the static-well model. 
We do this by modifying the standard Furry representation [2] through forcing the lowest- 
order Hamiltonian to be same as that used in that representation, but having the quarks 



in the proton provide the Coulomb field instead of assuming a point source. This requires 
the introduction of a new term in the interaction Hamiltonian we call the counter term, the 
effects of which however are quite simple to evaluate. We also define the static-well model 
and briefly review the calculation of the proton electromagnetic self energy. In section IV 
we use bound-state field theory to treat one photon exchange, and show that the results 
agree with the first order QM energy. We note here that in this paper we use the Coulomb 
gauge and treat only Coulomb photons for all exchanged photons. For the QM approach 
this corresponds to ignoring magnetic effects, and in field theory to leaving out transverse- 
photon exchange. In section V we then turn to two-photon exchange diagrams, which 
we break into two classes, one in which only one photon attaches to the lepton, with the 
other being emitted and reabsorbed in the proton, and a second in which each photon is 
exchanged between the lepton and a quark. In Section Va we treat the first class, which 
has no QM analog, and present a calculation of the contribution, which we call the proton 
Lamb shift. In Section Vb we treat the second class, but again make a breakup of the 
diagrams into firstly a part in which the proton is left unchanged, and secondly a part where 
it is excited. (In our model this means that in the spectral representation of the quark 
propagator it is either saturated with the Is state, or else that state is excluded). The 
first part will be shown to correspond exactly to the second-order QM energy. The second 
part, the proton polarizability, is then evaluated. The related calculation of the proton's 
static electric polarizability is carried out in Section VI, and it is shown that in the k, = 1 
angular momentum channel a complete cancellation between positive- and negative-energy 
state terms occurs, leaving only contributions from the k, = —2 intermediate states. In the 
conclusion we compare our results to the results of other calculations and describe directions 
for future progress. 



II. PERTURBATION THEORY OF FINITE-NUCLEAR-SIZE EFFECTS 

We consider a central potential for hydrogen or muonic hydrogen coming from a finite 
charge distribution p{x), normalized to unity. The corresponding static potential is 

V{x) = -Za [dx'j^^^. (3) 

\x — x'\ 



While we have Z = 1, the following discussion can also be applied to the case Z ^\. We 
start with a point-Coulomb binding field, so this distribution leads to the perturbation 

W(x) = -Za / dx'^-^ V^. (4) 

J \x — x'\ 

The first-order correction, 

E^^^ = f dx(Pl{x)5Vix)(P,{x), (5) 

is valid for either a relativistic or nonrelativistic calculation. We first consider the nonrela- 
tivistic limit. Then in leading order the wave functions may be replaced by their value at 
the origin and the first-order energy is 

E^'^ = -Za\MO)\' [dx fdx'P^^^l^^ 
J J \x-x'\ 

= -Za \M0)\'ldx'4n f-i-^x" + ..)j [p{x') - 6ix')] 

= ^ \U^)? j dxx'pix) = ^-^nvl j dxx'p{x\ (6) 

where in the last step we have assumed v to be an ns-state. 

The integral over dx has been carried out using a cutoff procedure which we now describe. 
We introduce a parameter u, understood to ultimately be taken to zero, and work with the 
basic identity 



/ 



-t „-u\x-Xl\ A 1 _ p-u\x2-XT,\ 

^^^^T r = ^V n' (7) 

\X2 — X\ \X — Xi\ W^ |iC2 — 33i| 

which for small u has the expansion 

1 e-"l=^~^il /I 1, , M, ,2 «^ ,s 

dx = 47r \X2 — cci H — \xo — xA \X2 — x-\\ + . 

\X2 — x\ |a; — a;i| \u 2 b 2A 

By differentiating once or twice with respect to u, we have 

Q-u\x-xi\ / 1 1 \ 

dx-. r = A7ci — --\x2-Xi\^ + 0{u)) (9) 

\X2 — x\ \U^ 6 J 

and 

Q-u\x-xi\ / 2 1 \ 

da; \x - xi\ = 4:71 { — - — \X2 - Xi\^ + 0{u) . (10) 

\X2 — x\ yw^ 12 J 

We have used Eq. ^ in the derivation of Eq. ([H]), and use Eq. flTUl) to evaluate the 
correction coming from the variation of the wave function to leading order. For S states, 
this arises from 

|0,(x)|2 = |0,(O)|2 (1 - 2Zam,\x\ + ...), (11) 



which yields an additional contribution of 



n{l) 



EY' = 2{ZaYmi |0, 



da; 1 33 1 / da:; 



, P{^') - Six') 



\X — X' 



2{Zafmi 10^ 



dx' A-K 



W 
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\x 



/|3 



+ . . .) [P(^') - Six')] 



27r{Za) 



■mi\(j)y{0)\'^ / dx\xfp{x) 



2{Za) 



mf / dxx^p{x). (12) 



This term cancels a corresponding term in second-order perturbation theory, which is given 
by the standard form 

(l)m[X2)(p in{Xlj 
f.. — 



^(2) = fdX2 fdXi(t)l{x2)6V{x2)J2 



6V{xi)(t)^{xi). 



(13) 



The sum over terms involving ?72 is — 1 times the reduced Green function. This expression 
is again valid for either a relativistic or nonrelativistic calculation. The correction to the 
potential is only non-zero outside the nucleus, which means that the wave functions and 
reduced Green function are evaluated for small arguments, because the Bohr radius for 
both the muon and electron is large compared to the nuclear size. We again take the 
nonrelativistic limit. Then the wave functions may be evaluated at the origin and the 
reduced Green function may be replaced by the nonrelativistic free Green function to give 



E, 



(2) 



'27r 



{Za)^mi 

2-11 



da:;2 / dx\ 



dx2 \ dx\ SV(x2) 

\U0)\' J dx', J dx', 
p{x'^) - 5{x'^) 



1 



\X2 - Xi\ 



6V{x^ 



1 



p{x',) - 6{x',) 



(14) 



3^2 — Xol *^2 — X\\ \X\ — Xi 

Introducing cutoffs allows us to carry out the integrals over the unprimed variables with the 
formulas given above, resulting in 

|2 



Ef = {Za)'mi\MO)\' fdx'.jdx'. 



X / da^i [^(ccg) — S{x'2)] \x'<2^ — xx 



PJO-Sjx',) 
\xi — x'-,\ 



(15) 



which yields 



E. 



(2) _ 7r(Za)^m; |0^ 



3n^ 



■m, 



/d../d..p(.,)|.,-.,|V.0-2/d.|.P,W 
/ da;2 / da:;ip(a:;2)|a;2 - Xi\^p{xi) - 2 / da; |a;|^p(a;) 



(16) 



As alluded to above, the second term in the square brackets in Eq. (fT6!) is cancelled by 
Eq. flT2l) . The first term in the square brackets is the third Zemach moment, which we 
denote {r^)z- An interesting feature about this term is that it too is cancelled by a term that 
arises when the nucleus is allowed to undergo low-energy excitations (proton polarizability), 
though we will not use this fact directly, and instead just evaluate the entire effect. We turn 
now to a field-theoretic approach based on the static-well model, and begin by introducing 
a formalism for bound-state field theory. 

III. FORMALISM 

While the formalism we use here is to our knowledge novel, it is a simple extension of 
the Furry representation [2], which we now briefly review. We will use this representation 
both for leptons and quarks, and begin by describing how it is used for the former. The full 
QED Hamiltonian used for describing the scattering of free leptons is H = Hq + Hi, with 
({^} = l,{c} = l,{e} = l), 

Ho= dx ip^x) [a-p + Pmi] iIj{x) (17) 

and 

Hi = q, J dx^{xh^^{x)A^{x), (18) 

with ge = ~e. We suppress normal ordering and the self-mass counter term for simplicity. 
The Furry representation is used when Hq is replaced by 



^0= f dxijU 



X] 



n Z^ 

OL ■ p + prrii — - — - 
\x\ 



ip{x). (19) 



This builds in a classical Coulomb field from an infinite mass proton. Carrying out a unitary 
transformation to eliminate Hq rather than Hq leads to the Furry representation in place 
of the interaction representation. While the interaction Hamiltonian Hi of a lepton with 
photons keeps the same form, the lowest order spectrum now consists of hydrogenic bound 
and scattering states, and the lepton Green function obeys the relation 

-ia ■ V^ + /3mi - -^ - z\ G{x, y; z) = 6{x - y), (20) 

which has the spectral representation 

Gix,y;z) = J2^^^^^^}^M, (21) 



In the present case we treat proton structure using the static-well model to specify the 
wave functions and Green functions of the constituent quarks. However, the influence of the 
proton on a lepton in a bound state cannot be treated perturbatively, so we need to build 
the binding of the lepton into the formalism nonperturbatively. We do this by modifying 
the breakup of the Hamiltonian given above to H = {Hq + Hx) + {Hi — Hx)- If we choose 

d.x^\x)-^^[x), (22) 

\x\ 

then this breakup is H = Ho + {Hi — Hx)- In the following we refer to Hx as the counter 
term, though of course it is to be distinguished from the electron mass counter term. We 
stress that we do not assume that a classical Coulomb field is present. However, because 
we have added a new term to the interaction Hamiltonian, Hq is unchanged from the usual 
Furry representation, and the same wave functions and Green functions used in Feynman 
diagram calculations in that representation can be used, although extra Feynman diagrams 
involving Hx need to be included. 

The use of the Furry representation must be extended to quarks in order to account for 
the proton's Coulomb field. In this case Hq and Hi are almost identical to the free case, but 
in Hq we assume the presence of a static well that confines the quarks, the details of which 
are given below, so that the up and down quark fields are expanded in terms of solutions to 
the Dirac equation in this well. The proton then consists of the usual two up quarks and 
one down quark, and the sum of their charges leads to the Coulomb field felt by the lepton. 

Since the lowest-order problem describes the basic physics of the atom, there has to be 
cancellation between diagrams in which a Coulomb photon is exchanged between the lepton 
and the quarks in the proton and diagrams with one counter term. In second order another 
cancellation between two-photon exchange diagrams and diagrams involving one and two 
counter terms must take place, and so on. Because the proton is now modeled as a finite 
object, the cancellation will not be complete, and we will identify the parts remaining after 
the near cancellation as the subject of this paper, proton structure effects. It is of course 
vital for this procedure to make sense that the cancellation not only takes place, but that 
the perturbation expansion converges. We will present results for the first and second terms 
in the expansion below. In determining the order of the expansion we note that the counter 
term Hx is of the same order as two His. 

The static-well model has well-known solutions which we show partly to establish nota- 



tion. We represent the solution to the free Dirac equation in a spherically symmetric well, 
centered at the same origin as used in Furry representation, by 

if2{r)x-K,,{r) 
Here Xkm(^) ^^ ^ spherical spinor and the radial wavefunctions obey the equations 

^ ^"'''/2(r) + [i?-m(r)]/i(r) = 0. (25) 



9r r 

Confinement is enforced by choosing m{r) to be constant inside the well and tending 
to infinity for r > R, which leads to the MIT bag model boundary conditions [v\. The 
ground-state solution for the case m{r) = for r < i?, which will be used throughout this 
paper, is 



fi{i") = — sin(wx) 

Mr) = - 
r 



sin[wx) 

coslwx) 

wx 



(26) 



with X = r/R and 



iV = 1 J^^ . (27) 

y/R\ 2w^ + cos2w; - 1 ^ ^ 

Here w = 2.042 787 for the ground state, which gives eg = 335.9 MeV for R = 1.2 fm. When 
this wave function is used for the up and down quarks the second and third moments, which 
can be calculated analytically, are 

, n, n 2W^ - 2W'^ +4:W -3 n , , 

and 

3 ^ 3 2W^ -2^7 + 3 ^ 3 

^ ' %w{w-l) ^ ' 

While an analytic form can be derived for the third Zemach moment, it is lengthy and 
involves the sine integral, so we give only its numerical value, 

(r3)^ = i?3 1.280 621. (30) 

The wave functions of the up and down quarks are identical in this zero mass case, and 
we denote them as 0g(ic). This approximation leads to the important simplification that 

10 



the charge density of the proton, even though it consists of three quarks, can be written in 
terms of 0g(a;), 

Pg(^) = 0g(a=)0g(a=)- (31) 

We will see that in our calculations of one- and two-Coulomb photon exchange this density 
will enter in exactly the same way as it does in Sec. II in parts of the calculation, thereby 
reproducing the results of that section with the static-well model charge density. Extra 
terms arising from the field theory approach will be identified with polarizability effects. 

Were we to use different masses, the charge densities of the up and down quarks would 
differ, and in that case we would use 

Pgix) = - (plix)(f)^{x) - - (l)\{x)(pdix). (32) 

The 2si/2 and 2^3/2 atomic wavefunctions for the lepton are the standard Dirac-Coulomb 
solutions and will be denoted as (f)v{x). As proton structure effects are strongly suppressed 
for the 2p3/2 case, even though we will continue to use v in formulas, in practice we will 
always assume v = 2si/2- 

The radius R is the only variable in this calculation, and we will use different values to 
study the R dependence of what is by far the largest proton structure correction, the finite 
size correction from one-photon exchange. However, because all other proton structure 
effects are much smaller, the value 1.2 fm is understood to be used for those corrections. 

For the calculation carried out here, which involves a quark propagating in the well, we 
use the same kind of spectral decomposition as given in Eq. ( 12T1) . 



G,(x.y;.) = ^ ^"'">^"''''> . (33) 

The sum over m for the lepton and n for the quark Green functions can be carried out using 
the method of finite basis sets 15|], which have been used extensively for atomic calculations, 
with only minor modifications of the associated computer code required for application to 
the quark Green function. This is because the atomic calculations were set up in the same 
kind of confining well as used here, but in that case only for the purpose of discretizing 
the spectrum, with the well radius chosen to be much larger than the atom or ion being 
considered. 

We will need the explicit form for the spin-up and spin-down proton wave function, with 

11 



the former being 



|Pt) 



^ijk 



-9h^ h^ h^ -I- h^ h^ h^ -I- h^ h^ h^ 



|0). 



(34) 



In Eq. ( I34l) a and h denote spin up and down states of an up quark and c and d spin up and 
down states of a down quark; eijk is the Levi-Civita symbol, which makes the proton a color 
singlet after the implicit sum over colors ijk is carried out. Because we are taking the up 
and down quark masses equal to zero, the associated wave function 0c(2;) can be replaced 
by (j)a{x) and (j)d{x) by ^^(a;), which simplifies later formulas. When the spin state of the 
quark is not important, we simply use 0g(a;). 

Energy shifts are calculated with the use of S'-matrix techniques, where we use Sucher's 
generalization of the Gell-Mann Low formula [16,], 

AE= lim -^\Yi(v\T\e-'^'''^'^\v). (35) 

Here Hi{t) indicates that a factor e~^'*' is included in the time integral over the Hamiltonian 
density in order to adiabatically turn off the interaction at large positive and negative times. 
The advantage of this formula is that the S'-matrix can be described with standard Feynman 
diagram techniques, with the adiabatic factors usually trivially leading to a factor 1/e that 
cancels the e in the numerator of the above formula, though when we deal with two-photon 
diagrams, the formalism is needed to cancel disconnected diagrams. Details of how this 
works along with other technical issues can be found in Ref. |3|. That work described a 
calculation of two-photon exchange diagrams contributing to energy shifts of excited states 
of heliumlike ions, but the basic approach is almost identical. The most important difference 
is that while in that work \v) describes two electrons, here it describes one electron or muon 
and three quarks, and is given by 

\v) = bl \p) . (36) 



The diagrams that involve one photon are given in Figs. [H [21 and 3(a) Fig. [T] is the 
standard one-loop Lamb shift, which has been evaluated to spectroscopic accuracy. Fig. [2] is 
the electromagnetic self energy of the proton, which was calculated using the Feynman gauge 



in Ref. 



^. The simplest diagram to evaluate is Fig. 2(a) , photon exchange between pairs of 
quarks. It contributes —0.222 MeV to the proton electromagnetic self energy, all of which 
comes from the vector part of the photon exchange. More difficult to calculate is the self- 



energy diagram of Fig. 2(b) This diagram in general requires the inclusion of a self- mass 



12 



counter term, though not for the zero mass case. After this subtraction an ultraviolet- 
divergent vertex term generally remains, but because the quarks move freely within the well 
the ultraviolet divergent part of this term vanishes, and numerical evaluation yields 0.658 

n 

MeV (the results given in Ref. |4l] referred to in the introduction are based on the radius 
R = 1 fm, but because they scale as 1/R a factor 1.2 must be inserted for comparison to the 
present work). We note that vacuum polarization terms do not contribute for an isolated 
proton, but may be of interest for muonic hydrogen, an issue that will be discussed further 
in the conclusion. Finally, Figs. 3a and 3b describe exchange of a photon between the lepton 
and the quarks in the proton and the first-order effect of Hx respectively, and we now turn 
to their numerical evaluation. 

IV. ONE-PHOTON EXCHANGE 



We evaluate Fig. 3(a) in the Coulomb gauge, and consider only Coulomb photon exchange. 



This leads to the energy shift 

AE,c = -a [ ^^ <pJ{x)u^^l{y)My)- (37) 

J i"^ y\ 

We have used our approximation of having the up and down quark wave functions being 
equal, so the sum of the contribution of the three quarks gives the single term 0g0g. A direct 
evaluation of this diagram for a state with principal quantum number n gives a result very 
close to —niia'^/n^, with the difference attributable to relativistic effects and the finite size 



of the proton built into our model. The associated counter term in Fig. 3(b) contributes 



AEx = a I ^(l)J{x)(t),{x). (38) 



For the 2si/2 state it has the value 



AEx = , , 39 

7^8(1 + 7) 



where 7 = yl — {Zay. 

Comparing the sum of AEic and AEx to E^^^ [Eq. ([5])], one sees that it is exactly 
reproduced by the field-theory expression. However, because we now have a specific model for 
p{x), we do not make the approximations made in the QM treatment and instead numerically 
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evaluate it. We find 

E(^) = 2.278 68 X 10"" a.u. [hydrogen] 

E(^) = 2.006 26 X 10"^ a.u. [muonic hydrogen], (40) 

where all digits are significant. The root-mean-square charge radius Vp for i? = 1.2 fm is 

rp = 0.728 966 i? 

= 0.874 760 fm. (41) 

We emphasize that this is not meant to be a prediction of the proton's rms radius, it 
is simply the static- well result when i? = 1.2 fm. However, if we use this in the standard 
nonrelativistic expression for the finite size effect in hydrogen given in Eq. (^ we get 

E^^^ = 2.211 Ih X 10"" a.u. [hydrogen] 

e'^'^ = 2.013 00 X 10"^ a.u. [muonic hydrogen], (42) 

which differ from the exact result by 0.07% and 0.34% respectively. Thus we see that 
the calculation reproduces the bulk of the nonrelativistic expression for the effect of finite 
nuclear size on 2si/2 energy levels. In the following, while the basic parameter of the static- 
well model is the well radius R, we use Eq. fHT]) to replace it by 1.3718 rp in all formulas 
dependent on R. 

As we have a complete model of the charge distribution, these small deviations can be 
attributed to higher moments and relativistic effects. By carrying out the calculation for a 
range of R around 1.2 fm, we find the fits 

^2i^(^p) = [2.97791 X 10""fj7 _5_i88Q9 x 10"i6f3]a.u. (43) 

for hydrogen and 

Eil\fp) = [2.631 70 X 10'^ fl^ - 1.126 29 X 10"^fJ|] a.u. (44) 

for muonic hydrogen, where fp denotes r^ in units of fermis (femtometers) . 

We first note that the coefficients of the first and second terms for hydrogen increase by 
close to a factor of {m^i^/me}^ and (m^/me)"^ respectively for muonic hydrogen, consistent 
with the dependence on mi shown in Eqns. ([6]) and (TT2|) . The coefficients agree at a level 
of the order of one tenth of a percent. We originally attempted the fit with a quadratic 
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term in fp instead of a term with the exponent 27, but were forced to use the latter form 
for hydrogen to get a proper fit. (The effect is less important for muonic hydrogen). In 
fact, it is known that the actual dependence of leading finite-size correction on fp is not 
the nonrelativistic quadratic form, but instead the relativistic form used above, as shown in 



Ref. [12|. In perturbation theory, the leading effect of the fractional power is a correction 
given by the Taylor expansion of r^'"^, which leads to a logarithmic term of relative order 

V. SCREENING OF THE PROTON SELF ENERGY 

As mentioned in the introduction, the proton electromagnetic self energy has contri- 
butions from the Feynman diagrams of Fig. 2. Before discussing how these diagrams are 
modified when the lepton interacts with the proton, which can be thought of as the Lamb 
shift of the proton, we mention another Lamb shift related term. This other effect, while 
negligible for hydrogen because it is of relative order {m^/M^y, makes a small contribution 
for muonic hydrogen, and is not suppressed at all for positronium, accounting for the self 
energy of the positron in that system. It was first derived by Fulton and Martin |13| , and is 
given by 



-EsEN - - — ^— 2" 
6Txn-^ m^ 



In I — ^ ) 610 -In koin,!) 



(45) 



This recoil effect, which shifts the 2si/2 energy in muonic hydrogen by 0.010 meV, is not 
included in our approach. 

An isolated proton can of course emit and reabsorb a photon, giving rise to the elec- 
tromagnetic self energy of the proton just mentioned. This contributes to the mass of the 
proton, but when the proton is in a bound state an additional shift arises, described in lowest 



order by the Feynman diagrams in Figs. 4(a) and 4(b) We note in passing that these dia- 



grams do not have an analog in the QM treatment, as they involve internal electromagnetic 
interactions in the proton. In this paper we restrict our attention to the second of these 
diagrams, which we refer to as exchange corrections to the EM self energy of the proton and 
label as Ai5ex(pLS). This is justified by the behavior of the lowest-order proton electromag- 
netic self energy, where the size of the exchange term and the quark self energy are of the 
same magnitude. The photon propagators are both taken to be Coulomb propagators, and 
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it is straightforward to show that this set of diagrams gives the energy shift 

loTT^ J \x — y\\w — z\ 

X J2 _ {p I • '^\i^)Mx)'4^^ jiy)My)<Piiz)'ipj (^) -Ip)- (46) 

In this equation i and j are understood to be summed over the two flavors and appropriate 
sums over color indices are imphcit; (ji is the charge of the up or down quark in units of e 
depending on whether z = 1 or 2 respectively; ipi, tpj, and their adjoints are field operators, 
but (pn and its adjoint are wavefunctions, with the sum over n going over all allowed values 
of angular momentum, /€„ and /i„, and the positive- and negative-energy states associated 
with Kn- For two Coulomb photons, one can show that only k^ = — 1 yields a non-zero 
contribution. After taking the normal-ordered product of the quark fields between the spin 
up wave function of the proton and noting the integrations over angles are all elementary, 
one is left with the sum 

AEex(pLS) = -a^^ / dxx'^Rggix) dyy'^Rgniv) 

[n] ^S ^" "^° -^^ 

/ dzz'^Rngiz) / dww^Ry^{w) -, (47) 

Jo Jo max{x,y) max{z,w) 

where [n] denotes summation over all s-states except the Is state, and 

R,,{x) ^ Mx)fi,ix) + f2i{x)h,{x). (48) 

We find, using finite basis sets to carry out the sum over [n], that 

A£;ex(pLS) = 0.5 X 10"^ meV. (49) 

The smallness of this effect is related to the fact the integral over z in Eq. fl47p is 

rR . 

/ dzz^R^^{z) -, (50) 

which vanishes when w > R from the orthogonality of the radial wave functions. This 
restricts the muon wave function to lie within the nucleus, which gives a large suppression 
factor that scales as {R/aoY- 

This same suppression factor should make the diagram of 4b, which is more difficult to 
evaluate, numerically unimportant. However, an interesting application of our approach 
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a 



would be the calculation of vacuum polarization. This was shown in Ref. j^ to vanish for 
a free proton, but when bound in an atom the arguments for its vanishing no longer apply, 
and in fact this should correspond to the effect of hadronic vacuum polarization, which is 
treated as a small effect, estimated in Ref. 5| to be 0.011 meV. However, we are not aware 
of any direct calculation of this term with zero-mass quarks, and will discuss how such a 
calculation could be carried out in the conclusions section. 

Before turning to two-photon effects, we give more details about the use of finite basis 
sets. As described in Ref. 15|], atomic finite basis set calculations are carried out in a 
well much larger than the atom, but with the same boundary conditions as used for the 
quarks. We continue to use an atomic basis set appropriate for hydrogen, but add a second 
basis set for the quarks, which is obtained by simple modifications of the atomic code, 
involving changing the fermion mass to zero, eliminating the potential, and changing from 
atomic units to MeV-fm units. Atomic grids are created on an exponential grid of the form 
'i"{i) = ^0 [e^'-*"^-' — l] . We use the same grid for both quark and lepton wavefunctions. This 
is done by choosing parameters such that if, for example, a 1000 point grid were used for 
the atom, the two hundredth point would be at r = 1.2 fm, so that the quark wave function 
would be put on a 200 point grid that matched the atomic grid, though of course the quark 
wave function vanishes for i > 200. Several grids were used to test numerical stability. As is 
also the case for leptons, a complete set of positive and negative energy states result for each 
possible value of k„, in this case A^ positive energy s-states and A^ negative energy s-states, 
with a typical value of A^ being 50. For leptons the effect of the negative-energy states is 
generally very small, entering at the order of the Lamb shift. However, for quarks they play 
a more important role. 



VI. TWO-PHOTON EXCHANGE 

We now turn our attention to diagrams shown in Fig. 5. In this section we will show that 
they in part reproduce the second order perturbation theory expression for E^'^\ Eq. (IT3l) in 
Section II, but have in addition extra terms we identify as proton polarizability. To compare 
with individual diagrams, it is useful to employ Eq. 1^ to represent 6V in terms oi p — 6, 
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which yields four terms for Eq. (fT3l) : 



i^l? = o? 



/d.di,4W±^ ^"'">^"'"> ±feW. (51) 

J \'^\ _L ^1} ^"^ \y\ 

/d.dj;dM,4W±^ ^"'">^'"faV"'V „W. (52) 



4? =-a' /dxdydz4(^)^E ^'"'"' ^"'^' ^^.(y). (53) 

and 

£g) = „^ /■ d. d„ d. d» ^l(.) -^ Y. ^-'"' ^"'^' r^ *.(»)■ (54) 

^'^ i \x-z\^ e^-em \y-w\ 

(2) 

The simplest diagram, Fig. 5a, is easily shown to give E\^ . The diagrams of Fig. 5b and 
its complex conjugate give E^^ and E^^ if we identify p{x) = Pg{x) as in the treatment of 

eaves c 
to it comes from Fig. 5c, given by 



(2) 

one-photon exchange. This leaves only Epp to be accounted for. In general one contribution 



AE2{ij) = a^qiqj^ 



dxdydzdw (j)l{x)(f)m{x)(j)l^{y)(f)^{y) 



x-z\\y-w\ e^-er. 



X (p |: ^\{z)i'i{z)i)]{w)i)j{w) : p^ , (55) 

where the definition of qi is the same as in Sec. El The various contributions to this term 
happen to sum to zero for the proton in our model, but we note they would not were we 
considering the neutron, or if we were taking the up and down quark wave functions to be 
different. 

The final two diagrams, 5d and 5e, are referred to as the ladder (L) and crossed ladder 
(XL) respectively. The closed loop in these diagrams is associated with an integration over 
a virtual energy. Because in this paper we consider only two-Coulomb photon exchange, the 
analysis of the loop integral is considerably simpler than the case where the photon prop- 
agators have energy dependence. That complication was encountered in the full Feynman 
gauge analysis of ladder and crossed ladder diagrams in excited states of heliumlike ions 
3j, upon which the present calculation is patterned. Issues involved in carrying out the full 
calculation will be discussed in the conclusion. In this simpler case we can carry out the 
integral over the timelike component of the loop momentum with Cauchy's theorem, which 
requires identifying the poles coming from the propagators. We partition the lepton and 
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quark propagators into positive and negative energy parts, which could lead to four contri- 
butions for each diagram, but the position of the poles is such that only two contributions 
fail to vanish. The surviving terms are 

A p /, ix ~2 2 [ dxdydzdw -^ 0t(a;)0^(a;)0]„(y)0^(7/) 

Z\£/l(++j = qi a / > ■ 

J \x-z\\y-w\ ^-^ e^ + e„-em-en 

{p\-i^lM)^n{z)i)^ai.'w)i^qr{'^)-\p) (56) 

and 



AEj^{—) = -q,Wj 



' "^ ' m-fi- ^ 



(p I : ^l {z)Mz)^iM^q. M : I p) (57) 

for the ladder diagram, and 

. , . .2 2 f dxdydzdw ^^ 0];(a;)0„,(a;)0l„(y)0^(y) 



X Uj Wy Z\ 6n- C^ ~r ^rn ^n 



{p\:iljl{z)-^n{z)-^i{w)iljgXw):\p) (58) 

and 



AExl{-+) = -qi^a 



~2„2 / da; dy dz dw ^ 0j;(a;)0^(a;)0|„(j/)0^(y) 

X Uj \\y Z\ 6g ^u ~r ^i^^ 6^^ 



(p I : ^J^ (2)^„ (z)^;^ (^)^g, {w):\p) (59) 

for the crossed ladder. The relative minus sign between the two ladder and the two crossed 
ladder terms comes from closing the contour in different ways, and we note that the Coulomb 
propagators are different for the ladder and crossed ladder diagrams. In the following we 
will give formulas only for AEi^{++) and Aii^xL( — I") as the other terms differ from these 
only by an overall minus sign and a reversed role of positive and negative energy states. 

A. Elastic contribution 

At this point we divide the sum over intermediate quark states into a part with n = g, 
which leaves the proton unchanged and corresponds to elastic scattering, and the remaining 
part, corresponding to inelastic scattering, which gives the polarizability contribution. We 
now show that the elastic term contains EpJ . As e^ is positive energy, only AEi^{++) 

19 



and AExl{ — h) contribute. Because the quark states are now all s-states, the angular- 
momentum dependence is particularly simple. The Coulomb propagators can be simplified 
with the replacements 





1 




-> 


1 


\x 


1 


z 


—)■ 


ria 

1 










y 


— 


w\ 




r2a 



(60) 

with ria = max(|a;|, \z\), r2a = iT^ax(|?/|, l"*^!)- (Pot later use we define rn, = min(|a;|, \z\), 
r2b = min(|i/|, |if |).) Integration over angles then yields 

J riar2a -^ e^-em 

with m being forced to be an s state. We have introduced the shorthand 



dV = dxxM dyy^ dz z' dww\ (62) 

Jo Jo Jo Jo 

Because i?gg(w)/47r = Pg{w), this reproduces the part of Epp in which the electron prop- 
agator involves sums over positive energy states. In order to include the part with negative 
energy states the crossed ladder must be considered. The minus sign mentioned above is 
needed for the energy denominators to match, and the fact that p{z)p{w) is symmetric 
under interchange of z and w is also needed. Thus part of the ladder and crossed ladder 
diagrams together with the other diagrams of Fig. 5 simply reproduce second-order pertur- 
bation theory with a particular form for the charge distribution of the proton. 

We evaluated these terms numerically as a test of the coding, since the QM treatment 
yields a nonrelativistic (NR) limit for comparison. The result for muonic hydrogen is 

E^^J = -3.948 X 10"^ a.u., (63) 

which is within 4 % of the NR formula. We recall that formula has a mixture of the Zemach 
term and a (r^) term, where the latter cancels part of the one-loop result. 

B. Proton polarizability 

We now turn to the evaluation of the remaining parts of the ladder and crossed ladder. 
The quark propagator is treated in the same manner as described in Sec. |Vl The feature 
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of the basis set mentioned in that section, whereby it automatically breaks into a positive 
and a negative energy part, makes separation of the various diagrams a simple matter. We 
begin with the case where n^n = i^n = ~^, but where we do not allow the intermediate 
quark to be in the ground state. The result for the ladder is 0.000 5 meV, and is almost 
completely cancelled by the crossed ladder, which contributes —0.000 5 meV, making this 
channel completely negligible. 

In order to treat higher partial waves we must now specify the spin state of the atom. 
Polarizability is important only for the 2si/2 state, which of course can be a triplet or singlet 
state once the spin of the proton is considered. For one- Coulomb-photon exchange these 
states have the same energy, as hyperfine splitting comes from transverse photon exchange, 
which is not treated here. However, with two-photon exchange, even when both are Coulomb 
photons, the triplet and singlet energies differ except for the / = partial wave discussed 
above. We therefore present in the following formulas for the energy levels of the 2si/2 
singlet and triplet. 

While each Coulomb propagator has its own partial wave expansion, leading in general 
to a double sum over /i and I2, the fact that the quarks are all in s-states leads to the 
simplification that these values are equal, li = I2 = I. Another simplification is that for a 
given /, only Km and /t„ values associated with that value give a non-vanishing contribution. 
We define angular factors Ai{Km, ^n) for the ladder and Bi^k^, ^n) for the crossed ladder 
diagrams, which are rational fractions resulting from the integrations over angles and sums 
over magnetic quantum numbers, tabulated in Table II for / = 1 and / = 2. Both the 
triplet and singlet state results are given, though in this paper we are only interested in fine 
structure. / = 1 corresponds to the dominant dipole transition. In terms of these coefficients, 
the formula for the ladder is 





[dv ""^^""^^ V 

la 2a m+n+ 


and for the crossed ladder 


is 




[dv ''''''^' V 

3a '4a m+n_ 



E4 / N Rvm{x)Rmv{y)Rgn{z)Rng{w) , . 
A{fim, l^n) (64) 



^v ~r ^ff ~r ^r; 



p , ^ Rvm{x)Rmv{y)Rgn{z)Rng{w) 

^v ^g ' ^n ^m 



where r-^a = max(|aj|, |iy|),r3b = min(|a3|, \w\) and r^a = niax(|y|, |z|),r4b = min(|2/|, \z\). 
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1 = 1 








1 = 2 


Km l^n ^l(T) Ai(S) Bi{T) Bi{S) 


l^m 


^n 


^2(T) 


A2iS) B2{T) B2{S) 


1 1 17/729 57/729 37/729 17/729 


2 


2 


26/1125 


66/1125 46/1125 26/1125 


1 -2 64/729 24/729 44/729 64/729 


2 


-3 


64/1125 


24/1125 44/1125 64/1125 


-2 1 64/729 24/729 44/729 64/729 


-3 


2 


64/1125 


24/1125 44/1125 64/1125 


-2 -2 98/729 138/729 118/729 98/729 


-3 


-3 


71/1125 


111/1125 91/1125 71/1125 



TABLE I. Angular coefficient factors for the allowed values of Km, Kn for / = 1 and / = 2. Factors 
of Ai are for the ladder diagram; factors of Bi are for the crossed ladder. "T" refers to the triplet 
atomic state, while "S" refers to the singlet. 

These can be evaluated with the techniques described above with only simple changes, as 
the spline basis set contains all values of k needed and the radial integrals over r[/r^+^ are 
evaluated with techniques valid for arbitrary /. The main numerical problem was ensuring 
that the typical rapid convergence of splines in atomic physics carried over to this problem. 
The calculation shows that the effect of the / = 2 channel is very small, with almost the 
entire effect of polarizability coming from the / = 1 dipole channel, which shifts the 2si/2 
energy down by 0.026 meV. This leads to our main polarizability result for the splitting. 



S(2si/2-2p3/2) = 0.026 meV, 



(66) 



consistent with the results from dispersion relation analyses, as will be discussed in the 
conclusions. 

VII. STATIC POLARIZABILITIES 

One of the basic electromagnetic properties of the proton is its static polarizability, which 
has an electric and a magnetic component. We concentrate on the static electric polariz- 
ability, which is given by the Particle Data Group as 



Op = 12.0(6) X 10"^ fml 



(67) 



While this value is extracted from Compton scattering data together with dispersion theory 
arguments [l8| , in the static- well model the lowest-order expression comes from the diagram 
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of Fig. 6, where the hne with a cross indicates a constant electric field with magnitude i?ext- 
The energy shift is related to the static electric polarizability through 

AE^ = -^apEl,. (68) 



We note two similar calculations using the closely related MIT bag model 19|, 20| have 
been presented, but detailed comparison with their results is not possible. A factor of Air 
has been absorbed into the external field. 

The basic equation for the energy shift of a proton in the presence of a constant electric 
field described by the potential 0(r) = — i^cxta^s is 

AEo. = qi^El,J2^- fdxdy{p\:ijlix)x3Mx)i^liy)y3i^gAy)-\p) (69) 

which reduces to 

AE = -E^ V /" d diy '^'^l(^)^30n(a;)0U^)y30a(y) + 2(t>l{x)x3(t>„{x)(t>l{y)y-My) 

(70) 



For this case the orientation of the quark spin does not matter, so the expression simplifies 

to 

"p = -2 ^ —^ — / da; dy 0j(a;)x30„(a;)0j^(7/)y30g(y), (71) 

„ eg en J 

where we have now isolated the factor —■^E'^^^ to give the formula in terms of Op. We 
emphasize at this point that this is a field theoretic derivation, and that the sum over n is 
complete, including negative-energy states. 
Introducing the notation [see Eq. fHHj) ] 

rR 



r. 



V 



dxx^Rij{x), (72) 



'0 

after integration over coordinate angles we find 

2_ ip ip A ^ ip rp 

X *- ' 0-n.' r)o- ^c X ^ ' err).' ■ 



_ __ Y^ gn' ng _ ^ V^ ' gn' ng /^^\ 



^ n ^g ^n y n ^g ^n 



Kn — — 2 



We calculated the two terms with basis set techniques. The most striking result found 
was that the k = 1 channel vanishes. If one sums over only positive energy states one obtains 
a nonzero result, but inclusion of the negative energy states leads to an exact cancellation. 
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This situation does not occur for the k = —2 channel, which in our model is solely responsible 
for the proton static polarizability. 

We also calculated a-g using a form for the quark propagator that does not rely on the 
spectral decomposition given in Eq. fl33|) . After rewriting Eq. f l7T|) as 



ap = 2 dxdy (j)l{x)x-iGq{x, y; eg)|/30g(2/), (74) 

we use the fact that the quark Green function satisfies the differential equation 

{-ia ■ V^ + /3mq - eg)Gq(a;, y; eg) = 5%x - y) (75) 

with appropriate boundary conditions. Without those conditions this is the same equation 
that a free fermion Green function satisfies, and a simple modification of the well-known 
partial-wave expansion of the latter propagator can be made to solve for the quark propa- 
gator. We illustrate this for the Green function of a massless quark with positive k = / and 
X > y/m which case 



X ( jKeg2/)x^^(y) -ijz~i(egy)x^-KM(i') ) ' ^^^^ 



The coefficients Ai are determined by the boundary condition. The fact that the propagator 
ranges over a finite volume allows the admixture of the solution proportional to Ai which is 
forbidden for a free propagator because of the boundary condition as x — )• oo. When /t = 1, 
the integrand in Eq. ([71]) includes the factor 

\-i]i{e^y)x[t,{y) I 

which vanishes. An analogous factor on the left-hand-side when x < y also vanishes. How- 
ever, for K, = —2 there is no such cancellation, and using 

we find the numerical result 

ftp = 25.4 X 10"^ fm^ . (79) 
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The spectral decomposition discussed previously gives the same result, and one can see 
the sum is dominated by the first P3/2 state, with higher-n positive-energy and negative- 
energy states entering at under one tenth of a percent. This result is a factor of two larger 
than the experimental value in Eq. fIBTl) . However, we note that the relation of the PDG 
result quoted to static polarizability determined from energy shifts involves some subtleties. 



discussed in [18] and more recently in [21|. 



VIII. CONCLUSIONS 

In this paper we present an approach to calculating the effect of the electromagnetic 
structure of the proton on energy levels of muonic hydrogen that uses a simple bound-state 
model for the proton, with only one free parameter, the radius of the well, which we have 
chosen to be i? = 1.2 fm. Once the formalism is set up, standard techniques of quantum 
field theory can be used to evaluate proton structure effects using this parameter, in contrast 
to the standard approach, which involves the analysis of forward photon proton scattering. 
We do not claim our approach is better, only that it introduces a different way of looking 
at the problem. 

We have considered only Coulomb photons in this paper. This is because our primary 
concern was in setting up the basic formalism, and testing it in the relatively tractable case 
of Coulomb photon exchange. To continue, while one could stay in the Coulomb gauge and 
introduce transverse photons, it is simpler to simply change to the Feynman gauge. Because 
the quarks all have the same energy, much of the work here is effectively already in the 
Feynman gauge, with 7^ ■ ■ ■ 7^ — )■ 70 ■ ■ ■ 7o- Continuing the calculation to this more complete 
state will allow the treatment of the interesting case of hyperfine splitting. The difficulties 
of proton structure are well known to be exacerbated in this case, with even ground-state 
hydrogen hfs uncertain at the sixth digit. Our approach should allow a direct calculation of 
what is often referred to as dynamic proton polarizability, as well as allowing a systematic 
treatment of other spin dependent effects that sometimes are difficult to disentangle. These 
calculations are complicated by the fact that the loop energy integral can no longer be 
evaluated with Cauchy's theorem. Instead one must carry out a Wick rotation to the 
imaginary axis and carry out the integral numerically. This rotation requires care because 
of the presence of poles and cuts in the complex plane, which leads to a number of extra 
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terms. 

We have shown that in our method a term, while extremely small, that can be thought 
of as the Lamb shift of the proton arises. The smallness of the effect had very much to do 
with the fact that all three quarks are taken to be in the Is state, which limits the k values 
allowed for the propagator. It is an interesting open question as to how this effect would 
change if corrections to the proton wave function involving non-s states were present, as 
they presumably are because of gluon exchange, but that is outside the scope of this paper. 

We have shown that the one-photon diagrams together with the part of the two-photon 
diagrams where the proton is unchanged give the result 

Efns (non-recoil) = [-7.161 22 f^^ + 0.046 70 fj] meV. (80) 

This was calculated without the large recoil corrections present in muonic hydrogen: mul- 
tiplying the appropriate factors of {mj./m^)^ and {rrij. / m^j,)'^ into the first and second terms, 
our result becomes 

Eins = [-5. 199 5 f^T + 0.030 5 f^] meV. (81) 

The polarizability of the proton was calculated in a novel fashion. The basic result of this 
paper is that despite the difference in approach, a similarly small result is found. Specifically, 
again restoring recoil corrections by multiplying our result by (?7ii./m^)^, we have found the 
proton polarizability correction 

AEpoi = 0.017 meV. (82) 

While inclusion of transverse photons may lead to quantitative changes, we consider it 
unlikely that a qualitative change will result. 

We now compare our results with Ref. [51], where the formula 

AEls = 206.0336(15) - 5.2275(10)rJ + AEtpe (83) 

is given in their Eq. (32), with the breakdown of Aii^xPE = 0.0351(20) meV given in Table 
I along with results from other calculations that differ by under 2 |ieV. 

The quadratic term compares well when a term —0.0275 arising from radiative corrections 
to finite size, not treated here, is removed. The breakup of Aii^TPE involves three terms, 
a subtraction term of —0.002 meV, an elastic term of 0.023 meV, and an inelastic term of 
0.014 nieV. The first has no counterpart in our calculation, but if we identify the elastic 
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term with our Zemach contribution of 0.020nieV and the inelastic with our polarization of 
0.017 meV we see fairly good agreement. 

Perhaps the most interesting calculation left undone in this framework is vacuum polariza- 
tion. This of course dominates the muonic hydrogen 2si/2 — 2^3/2 splitting when an electron 
is in the vacuum polarization loop, but in our framework we could also put in the zero-mass 
quarks confined in the well into the loop. For an isolated proton, symmetry arguments ^ 
show that vacuum polarization vanishes, but once in an atom the arguments no longer hold, 
and a finite effect should be present. The standard approach is to introduce pion loops, in 
which case the previously mentioned very small "hadronic vacuum polarization" value of 
0.011 meV results, but this approach is quite different. Such calculations are, however, par- 
ticularly challenging because of the high degree of divergence present, which always presents 
difficulties for bound-state methods. We are investigating whether techniques that have 



proved useful in studying vacuum polarization effects in atoms 22| , where careful grouping 
of angular momentum contributions allows an accurate treatment of highly- divergent terms, 
can be extended to this novel case. 
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(a) Self-energy (b) Vacuum polarization 

FIG. 1. The two diagrams contributing to the Lamb shift in the hydrogenic bound state. 



(a) Exchange (b) Self-energy 

FIG. 2. The two diagrams contributing to the Lamb shift in the proton bound state. 
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(a) One photon exchange between the lepton (b) Counterterm contribution 

and a quark 

FIG. 3. Order a contributions 
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(a) Proton Lamb shift vertex correction (b) Proton Lamb shift exchange correction 

FIG. 4. Order a^ corrections to proton Lamb shift 
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(a) Two counterterm diagram 



(b) One counter term one photon exchange 
diagram 
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(c) Zero-loop two photon exchange term 



(d) Ladder diagram 
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(e) Crossed ladder diagram 



FIG. 5. Order a^ contributions to polarizability 
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FIG. 6. Diagram for static polarizability 
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